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■ Abstract 

An influence of the weak microlensing effect on the pulsar timing is investigated for pulsar 
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B1937+21. Average residuals of Time of Arrival (TOA) due to the effect would be as large 
as 10 ns in 20 years observation span. These residuals can be much greater (up to 1 ms in 20 
years span) if pulsar is located in globular cluster (or behind it). 
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\q ! 1 Introduction 

O . 

First, the problem of the influence of the weak microlensing effect on the pulsar timing observations 
was discussed in (Sazhin, 1986). It was considered as interstellar Shapiro effect. The massive body 
that flies not far from the line pulsar-observer produces changes in the observing frequency of 
Qh the pulsar similar to glitches. Estimations were made for Crab and Vela pulsars, glitches in 

these pulsars can be partially explained by the influence of the effect. Substantial contribution to 
the problem was made by ( |LarchenkovafcDoroshenko, 1995| ); they mainly investigated the case of 
microlensing (i.e. the gravitational deflector flied very close to the line observer-pulsar). It was 
shown that the microlensing effect would cause short-term growth of the residuals and follow-up 
relaxation. Whole interaction would take less then several years and the maximum amplitude of 
the residuals would be 20-30 ms. Such remarkable events are very rare, but all the pulsars are 
affected by the weak microlensing effect to a greater or lesser extent. This effect was considered in 
( |Ohnishi et al.,1996| ), where timing of millisecond pulsars was proposed as detection method for 
MACHOs. Growth of number of observed pulsars and time span of observation would make such 
detection easier. Numerical estimates were made in | |Hosokawa, Ohnishi, Fukushima, 1999| ). They 
stated that even when the measurement accuracy reaches to 10 ns, probability of the remarkable 
influence would be in the order of 10 _1 for the pulsar of a few kpc distance from us observed over 
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ten years. On the other hand there's well developed formalism for the effect that came from the 
optics. The weak microlensing effect causes distant sources like quasars from ICRF to "tremble" 
on the level of tens of mas. It was shown in ( |Sazhin, 1996| |Sazhin, Zharov, Kalinina, 1998| ) that 
these angular fluctuations range from a few up to hundreds of microarcseconds and this leads to 
a small rotation of the celestial reference frame. In (Sazhin, Zharov, Kalinina, 2001) influence of 
the effect on parallax measurements was considered-apparent parallax can be even negative due 
to the influence of the effect. Also, the weak microlensing effect can affect VLBI observations 
(Sazhin, Pshi rkov, 2005| ) and it should be taken into account with new generation of space-based 
VLBI. In dKalinina, Pshirkov et aL~ 2006) some statistical studies with toy-models were made, that 
was applied later to real model of the Galaxy. In fact, both weak microlensing effect and fly- by 
effect on timing are very similar and can be considered as manifestation of 4D (four-dimensional) 
astrometry piyasov et al., 1 990) 

In this work we tried to apply eikonal formalism that was developed earlier for investigation of 
weak microlensing effect for use in pulsar timing studies. 

The paper is organized as follows. In Section 2 we give a short review of influence of a passing 
body on pulsar timing in eikonal approximation. In Section 3 we apply a model of distribution of 
stars in our Galaxy to numerical estimations of their influence on pulsar timing and conclude our 
consideration in Section 4. 



2 Eikonal formalism applied to pulsar timing 

Change of phase during the propagation of electromagnetic wave can be obtained as a solution of 
Hamilton- Jacobi equation for a massless particle: 

9 -|^ = o (1) 

y dx^ dx" v ' 

Though S formally is a function of action, we hereafter identify it as eikonal or wave phase along 
the trajectory of the ray of light. 

In weak field approximation the metric tensor of gravitational field can be written down in a 
following form 

9fj,u = r\^v + h^v (2) 

Here 7/ M „ - is flat Minkowskian metric, h^ v - small additions to the flat metric that describes 
gravitational field of spherically symmetric body (star) 

Equation (JJJ can be solved in the following form: we take an exact solution ( Weinberg, 1972| ) 
and then take its asymptotic when the impact parameter of the propagating ray is much larger 
then the Schwarzschild radius r g = J; 1 (M-mass of the deflector) 

i// = i/,l+ r -^-arch{-) (3) 
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Here ip is full change of the phase along the trajectory, ipi- change of the phase along the trajectory 
that corresponds to the propagation in the flat space and time, r g - Schwarzschild radius of the 
deflector, cj-frcquency of the electromagnetic wave, r- some point on the trajectory, p - impact 
parameter (i.e. minimal distance between deflector D and curve of photon propagation). 




Figure 1: Schematic picture of origin of the weak microlensing effect. Z?-dcflccting object. 

Only the second term in is a matter of interest to us, though it's only a small addition 
to the usual change of phase during the propagation. The complete phase shift can be obtained 
as a sum of two solutions. The first is a phase shift during propagation from the source of the 
electromagnetic waves (which is located in Ti n ) to the closest approach to the deflector (we set 
the point of origin to the center of the deflector): 

5ip- = i/j(r = r in ) - ip(r = p) = ^-arch( — ) (4) 

c p 

The second - is a phase shift during propagation from the closest approach to the deflector to the 
observers (at r out ): 

Sip- = ip(r = r out ) ~ ip(r = p) = ^^arcft.(^£) (5) 

c p 

And the total phase shift is : 



3 



Aib = -^—ln( - — ) (6) 

c p z 

We treated the deflector as a motionless body in this solution. In fact, all stars, including 
MACHOs of our Galaxy are moving. Approximate solution of space-time metric in the case of 
moving deflector and trajectory of photon in such a variable gravitational field was calculated in 
| |Kopeikin, Schaffer, 1999| ). The metric which originates from a moving body and small perturba- 
tions of photon trajectory in gravitational field of this body, differ in (-^) 3 terms from our solutions 
and we will omit this difference. To describe the motion we take p (impact parameter) as function 
of t only: 
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Indices denotes values at different epochs t\ and t%. The first two terms are negligibly small, 
so we can rewrite expression (JJJ: 



A^-A^ 1 = ^ln(4) (8) 
c P 2 

Also we can write out time dependence of pit): 
pit) = y/ Pq + v 2 it — to) 2 , here po-mmimal impact parameter, v- velocity of relative motion of 
pulsar and deflector, to - epoch of the closest approach. We can rewrite the equation for the phase 
shift and obtain equation for time delays or residuals of Time of Arrival (TOA). It's worth noting, 
that these delays don't depend on frequency of electromagnetic wave: 



ST = ^ ~ A ^ = r JL ln(4) = -^H i +V l? 2 ~ t0 t ) 



We can set the first epoch t\ equal to zero and discard the second index, t% = t : 

5T=- r JL H ti+f} t -f ) (9) 

Here, t is time span of observations (we set the epoch of initial observations equal to 0), to - is the 
epoch of the closest approach of the deflector to the line of propagation. 

It's convenient to consider this problem on the "plane of deflector". Thus we convert all linear 
measures into angular ones: po — ^od , v — l^d, 8q -angular distance of the closest approach of 
deflector to pulsar, p - angular velocity of the relative motion (mainly due to the proper motion 
of pulsar), d- distance between the deflector and the observer. Hereafter phrases like "deflector s 
close to pulsar" mean we observe close angular coincidence of the bodies, not in 3D space. 
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Apparent pulsar's position on the 
deflector plane 



Figure 2: The plane of the deflector: red circle represents apparent position of the pulsar. Deflector 
with proper motion /j, and impact parameter 9q is passing by near this position 

Value 8q depends on location of pulsar in Galaxy and its proper motion. The higher is density 
of deflectors in the neighborhood of pulsar on the celestial sphere, the smaller that value would 
be. We take into consideration only deflectors between the pulsar and the observer, because they 
make the largest contribution on the effect. 



We chose two pulsars J1643-1224 and B1937+21 for further estimates, because they're quite distant 
and located in populated regions of our Galaxy (B1937+21: G; = 57.51Gb = — 0.29r p = 3.6fcpc; 
J1643-1224: G/ = 5.67G fc = 21.22r p = 4.86fcpc) ( |ATNF-psrcat, 1999| [Manchester et al., 2005| ) , 
so probability that effect would have place is much higher than for other millisecond pulsars. It's 
essential to define values 9q and t e - average duration of influence. They can be approximately 
found in such way ( |Kalinina, Pshirkov et al., 2006| ): stars are nearly uniformly distributed in the 
neighborhood of the pulsar on the celestial sphere; the angular distance to the nearest star, which 
would affect the pulsar timing depends on the location of pulsar. 

We calculated the density of the stars in the neighborhood, using accepted model of the disk 
of our Galaxy ( |Bahcall, 1986| ). 



N(8, </>) - sought density in the direction of the pulsar, which is assigned by the angles #,</>. 9- 



3 Estimates for B 1937-21 




(ii) 
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angle between the line of sight and the Galactic plane, <f>- angle between the projection of the line 
observer-pulsar to the galactic plane and the line Solar system-Galactic center;^- distance from 
the observer. 

n(r,z) = n exp(-^^)exp(-^) P c- 3 

no = 0.1 -density of the stars in Sun's neighborhood, r -distance from the axis of the Galaxy, 
z-distance from the Galactic plane , Rq — 8000pc -distance between the Solar system and the 
Galactic center, 3500pc and 325pc - radial and vertical scales of the model, accordingly. 

r(f , 9, <j>) = \Jnl+(, 2 cos 2 (<9) - 2i? f cos(0) cos(0) 
Z (£,M)=£sin(0) 



JV(0,0) = /n o £ 2 exp(f) 



exp(- 



2 {B)-2R (, cos(0) cos(0) 



(12) 



3500 



)exp 



■ gsin(g) 
325 



R) 



Average angular distance 9\ between the pulsar and the closest deflector (star) can be found 
with taking into account N(0, <f>) . Values 9o and t e were calculated using Monte-Carlo simulation: 
a circle of 9\ were circumscribed around the pulsar on the celestial sphere, then a large amount 
(1000) of test deflectors with proper motion u were started from this circle under random angles 
a. As a result we found distributions for values 9q and t e , and their averages, that were used in 
following estimates. Only known distribution of stars in our Galaxy was used in our estimates and 
if we take into account possible influence of Dark Matter, then sought values can be lower in 2-3 
times, because mass of DM doesn't exceed mass of ordinary matter more than 4-5 times. Also, we 
set mass of deflectors equal to Mq. Values that are essential for further estimations (J1643-1224, 
B1937+21) are given in the table below. 



PSR 


9i 


M 


9o 


t e 


J1643-1224 


7.3" 


~ lOmas/yr 


4.7" 


470 yr 


B1937+21 


2.5" 


~ lOmas/yr 


1.5" 


150 yr 



We can see the influence of the effect on the residuals, but only trends of cubic order and higher 
will survive during usual fitting procedure (Bakcr&Hellings, 1986). Linear and quadratic terms 
will redefine apparent period of pulsar P and its first derivative P and can't be found. 

Residuals of TOA due to the weak microlensing effect can be written as follows: 



ST, 



post fit 



= Ct 3 + Dt A + Et 5 



(13) 



C, D, E are coefficients in Taylor's series of function (5) where t = 0. 

Plotted coefficient C depending on to is represented in fig. 0] (plotted for B1937+21). 
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One can see from the plot that the fastest increase of residuals takes place when the epoch 
of the initial observation are 50-150 years away from the epoch to, because the third derivative 
have maximum in that interval maximal. If the initial observation coincides with the closest 
approach, then only fourth and higher orders term will affect timing and the residuals will increase 
much slowly. Magnitude of the residuals after subtraction of linear and quadratic terms can be 

" +M 2 (t-t c " 



expressed as follows: ST postfu = - r -f ln( 9p+ ^ J o) ) - A(0)t - B(0)t 2 , where ,4(0), 5(0) -linear 



and quadratic coefficients at t = 0. 



A(0) 



B(0) 



_ 2r g jjL 2 t 

c e 2 + nHi 

, 2r g nH 2 



c 9 2 + c (6» 2 + y?tl) 2 

The plot in fig. shows magnitude of the residuals at different to (0, 50, 100 years; blue, green 
and red graphs accordingly). Module of that magnitude depends only on module to. 

Residuals of 10 ns magnitude due to the effect of weak microlensing will appear with probability 
of ~ 50% if time span of observations exceeds 20 years. 

We can also calculate Allan variance (AVAR) for pulsar time scale with time residuals caused 
by the effect. 

TO A residuals due to the effect can be significant, if Q\ (angular distance between the pulsar and 
the nearest affecting body) is much smaller than average. The plot Fig. 7 represents situation when 
9o = O.lmas. This situation has ~ 0.5% chance of probability in case of B1937+21; Probability 
reduces like 6q inverse squared. The magnitude of the residuals can be as a great as 800-1000 ns 
in the same 20 years span. 

However, if we used in fitting procedure terms of cubic and higher orders, then the magnitude 
of the effect can be effectively set to 0. 

The magnitude can be much greater for pulsars in GC (or pulsars behind GC) (n = 10 3-4 pc -3 , 
L (length of path of ray in GC)= lQpc, <i(distance to GC) = 1 — lOkpc)). t e and &o can be much 
smaller because the density of stars in GC is large, The magnitude of the effect will be much 
greater ( the same 1 ms in 20 years span) . Time of one significant interaction will be quite small 
(20-30 years). Complete investigation of the question can be found in ( |Sazhin, Saphonova, 1993, 
|Larchenkovafe Kopcik m72006| ) . 



4 Conclusions 

So, we can make several conclusions: average TOA residuals due to a weak microlensing effect is 
about 10 ns (B1937+21) in 20 years span. TOA residuals can be effectively set to zero by using 
higher order terms in fitting procedure (not for pulsars in globular clusters. Residuals can be 
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much greater if pulsar is located in a globular cluster, so the pulsars in globular clusters can't be 
recommended for using in PT scale. 
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Figure 3: TOA residuals, caused by the effect; no fitting conducted yet. 
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Figure 4: Graph for the coefficient C(t) shows that the influence of higher-power order items 
should be taken into consideration. 
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Figure 5: TOA residuals due to weak microlensing effect. The blue curve corresponds to to =0, 
green one to to=50 years and red one to io=100 years. 
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Figure 6: Allan Variance (AVAR) which appears from a weak microlensing effect. 
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Figure 7: Strong influence of the effect. TOA residuals can reach 1 ms in 20 years observations 
span. 
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